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Design Considerations for Radiometer Based on Uncooled 

Microbolometer Detector 

Ernest Grimberg -OPGAL 



Abstract 

This paper concludes the large effort st^ed by OK>a. in f '^j^l^jS^'lSXSri^^^^ 
is not a tadiometric one. by using the NUC flag as an extremely low frequency chopper. 

Keywords : FLIR, Temperature measurement, UDCOoled , Microbolometer .Radiometer. Real time. 



1. Introduction 

1.1 



Ooeal manufectures a large number of uncooled microbolometer FLIR cameras covering a large number of 
Xt^^^clerJ contain regular uncooled microbolometer detectors ^t do «>-J^7,-^;"j;^,^^^^ 
liM A laree amount of these cameras comes with a standard lens of 35mm or 50 mm focal lengtfi. The mam idea 
f to uptS^^r™ L order to provide radiometric temperature measurement capability under the restncUon 
of very minor hardware modifications. 
1 7 The tjaner presents the studies performed for developing the measurement concept. The paper presents a relatively 
]^ed mJS ?S Z> energy exchanged among the microbolometer detector and the view, the optics, and^e 
ttill cXia ZTa mSmatical model that describes the detector signal versus target temperature was deyi^. 

»n,TSuK a Mhmomial model for tempei^re measurement The paper also presents die solution 
mZ^L^^eLLg'Srm^ drift and the accuracy obtained. There is one mapr problem n« 

Ti^^Tm^ vtSr S two-dimensional modulation transfer fimction (MTF) influence on the temperaftire 
rei« hLTbe cU^^^^ by an inverse fUter. THis subject is only mentioned in this paper but . not 

deeply treated. 



2, Energy exchanged model 

2.1 The video signal at the microbolometer output is generated by the energy exchanged between the detector's 

elements and: 

The view 
• The optics 

\e si«nSS^u5?^lpTaTansient process. Assuming that the readout process is intemipted. the energy exchanged 
^ts pCS AeSy state situation' is achieved. The energy exchanged process time constant is equal to the rano 
between the suspended structure Thermal Capacity and Thermal Conductance. 



2 2 Assuming an unlimited time interval and assuming tbat the view can be modeled as a Blackbody radiation source 
(eniissiviQr=l) the energy exchanged between the detector and the view is expressed by the following formula: 

Equation 1. Energy in joules, exchanged between the detector and the view. 
P _ view represents the energy [in joules] exchanged between one detector element and the scenery. 
1(1, TempQ + ATemp) represents the view, modeled as an average Blackbody radiation source in Watts/(srra' ). 
TempQ represents the detector's suspended structure temperature in degrees Kelvin. 
A det represents the effective area of a single detector's element in m . 
/ # represents optics f number. 
A represents the wavelength in meters. 
A J represents the lower integral limit. 

^J^ represents the upper integral limit. 

Trans _optics{X) represents optics transmittance 

£{X) represents the emissivity of detector's suspended structure. 

t presents the time in seconds. 

Solving Equation 1 for the time variable, ±e energy exchanged between one detector's element and the view is described 
by the following expression: 

P vie^^Temdi=^^^^-^^^'^^^ -^ATemp)-IiAJemp,)y Trans _opHcKm^y 

Equation 2. Energy in joules, exchanged between the detector and the view. 
Figure 1 describes Ae energy exchanged between the detector and the view assuming: 

• f number equals 1 

• high transmittance optics (93%) 

• detector emissivity equals 0.8 

• effective detector area equaling (41 *10'^^ meters 

• a spectral band ranging from 8*10"^ to 12.5*10"^ meters 

2 3The energy P^optics(ATemp) exchanged between the detector and optics assuming simUar conditions to those described 
in Paragraphs 2. 1 and 2.2 is given by the following mathematical expression: 

P_opticK^emp)=nT^A_dcHsm/^^^^^ 

k_\ 

Equation 3. Energy in Joules, exchanged between the detector and optics. 
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Figure 1. Energy exchanged between detector and view. 

LiXJemp, + ATemp) represents the optics, modeled as an average Blackbody radiation source in Watts/(sr.m^ ) and 
emissivity equals (1 - Trans _ opticsiX)) . 

0 = flrc_tan(Y^) 

0 represents half of the planar cone angle related to the light collection. 

0 = 0-464 radian for 1 . 
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Figure 2. Energy exchanged between the detector and optics. 



2 4 The energy P camera(ATemp) exchanged between the detector and the internal camera parts assunung similar 
conditions to Aose described ta the previous paragraphs 2.1, 2.2, and 2.3 is given by the followmg mathematical 



expression: 



X 2 



/>_c«wer<Ar^m;7) = ;rr.>4_det.(l--(sin(^)') jmJemfi, ^LTemp)-Ii^Jemp,)) s{X)dX 



LiAJemp^ + Afemp) represents the internal camera parts, modeled as an average Blackbody radiation source in 

Watts/(sr-m^) and emissivity equals one. 
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Figure 3. Energy exchanged between detector and internal camera parts. 

The energy exchanged between the detector and internal camera parts is the largest one, followed by the 
exchangS^between the detector and the optics. The energy exchanged between the detector and the «ew b the smafl^t 
one. Assuming, shnultancously, the same temperature difference between the detector and *e;.ew, defector and op^^^^ 
md detector iainst camera internal parts, the amoimt of energy exchanged between the detector and the view ,s about 

2 5SS^?r,^s!I!?obtained shows that that the capabiUty to reconstruct the view temperature from the vWe^giuJ 
■ collated is very difficult. The energy exchanged by the detector with the ^^^^^"f^^^^'^^t^;^^^ 
SLtribution of fee internal camera parts. This distribution is not constant, due to the feet that the camera s Aeiml t^ne 
corwtant is very long (about half an hour) and usually the environmental conditions are not stable for such long penods 
of time. A dynamic model based on the temperature distribution of internal camera parts is too <=°n»Plf ^^f.^/y*"^^ 
Sie c^pabiurto be implemented in the existing hardware. Therefore a different solution has to be adopted. The mam 
Sea Tto ul ?he non'uniformity coixection (NUC) update pnK:ess in order to eliminate the camera paj 

influence. The flag inserted in front of the detector during the NUC update process is used as an extremely low- 
Squen^ ch^per The flag's surface facing the detector was treated and coated in order to behave as close as possible 
a? i S Ebody surfece. A few vely accurate thermistors {±0.025X) have been installed at the followmg 

locations: 

2.5. 1 One thermistor on the back side of the flag surface 

2.5.2 One thermistor on the external surfece of detector's vacuum package 

2.5.3 Two thermistors have been glued to the external optics metal case. 



Energy exchanged idtio View over Total 



0J24 



Pjview(Arernp) 



{P_camera(ATCTra))+P_qptics{ATanpHP^ 0.122 



0.12 



0.118 



























/ — 









0 



5 10 

Ternpctdtuoro iMerencd 



15 



20 



Figure 4. View energy exchanged versus total energy exchanged. 



3. Temperature measurement description 

S.lLet us assume, just for simplicity, that the optics influence is negligible. (Later on we shall review this assumption.) 
Since the NUC update process is periodically performed, we assume that during Process the average detector signal 
Asn (after hRJC update and bad pixel replacement (BPR) have been performed) has been calculated. We also assume that 
during the same NUC updating Process number yk the temperature of all the installed thermistors have been acquired. 
Tfn represents the temperature measured on the flag during the n NUC update process , 

?V„ represents the temperature measured on the external surface of detector's vacuum package during the number n NUC 
update process, 

T On represents the average temperature measured on the external optics case during the number n NUC update process. 

Let us assume that the average temperature of some arbitrary, relatively large area of interest (the MTF problem* will be 
mentioned later on) has to be measured. Let us fiirther assume that the specified area belongs to an ideal Blackbody 
radiation source. 

The first approximation of the measured temperature is expressed by Taylor series: 

T_mes=T^me^As)+ \ Sig-As)^~ -^^^ \ sig~Asf ^^^^ ^^^^ \ Sig-Asf 

were: T _ mes represents the first approximation of the measured temperature. 

Sig represents the average video signal after NUC and BPR of the area of interest mentioned above, at 
some point in time between the n and NUC update processes. 

We define As to be the average detector signal value measured on the flag during the number n NUC updating process. 
Therefore by definition: 



T _mes(,As„) = Tf„ 
T_mes = FiSig-As„)+Tf„ F(Sig-As„) = f^a, -^Sig-Asj' 



Equation 4. First approximation of the measured temperature 

where: 

F represents the transform that translates the video signal to temperature. 
Equation 4 presents the general idea of using any chopper as a reference source. The chopper*s video signal is subtracted 
from the video signal to be measured. In this way a differential measurement is performed. The video value 

Sig - As„ is translated to temperature and the offset value Tf„ is finally added. This general approach works very well 

in the case that As„dxid Tf^art almost constant between successive samples. However, as explained above in 

paragraph 2.5 this is not the case. The Error between the real object temperature and the reconstructed one by the first 
approximation formula, was plotted versus the temperature measured on the external surface of detector's vacuum 
package (Td) during a period of more than 10 minutes. See Figure 5. The plot described in Figure 5 shows a typical 
behavior. Figure 6 describes the minimum mean square error linear fit to the Error value. However, the slope of the 
linear best fit is not constant in time. A high slope value means that the camera is far firom thermal steady state situation 
and a very low slope means that the camera is in thermal equilibrium. 

3.2 Each time the NUC update procedure is performed, the error of the measured temperature Tjnes caused by the 

internal camera temperature fluctuations is eliminated. The typical behavior described in Figure 5 and Figure 6 suggests 
that by adding an additional term to the first approximation ( Equation 4) a better solution will be obtained. 

AT* 

T_mes = F(Sig-As„) + Tf„+ - " {Td - Td„ ) 



Equation 5 

where Tdn represents the temperature measured on the external surface of detector's vacuum package during the n 
NUC update process. 
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Figure 5. Temperature error versus vacuum package thermistor temperature. 
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Figure 6. Linear best fit to the Error value. 



The general idea expressed in Equation 5 cannot be directly implemented. An approximation of this idea was 
implemented using the values i45„ . % and Td^ producing Equation 6: 



Td„-Td„^, 

Equation 6. Second Approximation. 



The expression FiAs„ - As„^^ ) estimates the total flag change in temperature measured by the detector. However this 

value contains the error caused by the internal camera parts change in temperature and the real change of flag's 
temperature. 

3.3 This paragraph deals with F function approximation, A mathematical model that describes the detector's signal after 
NUC and BPR processes was built, in order to estimate the capability to translate the signal to temperature. Figure 7 
describes the output of the mathematical model. This mathematical model was used in order to investigate the accuracy 
that can be obtained using minimum mean square error polynomial approximation to the modeled video signal. For a 
relative small dynamic range of about 60 degrees Celsius a second order polynomial expansion gives an acceptable 
accuracy. For example, a second order polynomial expansion is adequate in order to model the optics contribution to the 
video signal. The third order polynomial approximation covers, with an acceptable error, a range of 300 degrees Celsius, 
while the fourth polynomial approxunation spans a range of about 400 degrees Celsius. 




Figure 7. Video output mathematical model 



The minimuin mean square error polynomial expansion of order L that samples the data at Q different values, requires 
the solution of i +1 Hnear set of equations. However, the arithmetic used in order to solve the linear set of equations has 
to be carefully chosen. The minimum number included in the linear set of equations is Q , while the largest number is 

Q 

Sig^ . The graphs presented in Figure 8» Figure 9, and Figure 10 have been calculated using high arithmetical 

1=1 

accuracy. A practical solution might be to divide the temperature range into a number of small regions (including some 
small overlap between the regions), using a lower-order polynomial expansion at each region. Up until this point we 
neglect the optics temperature influence. The next paragraph describes the solution adopted by Opgal in order to deal 
with the optics temperature influence. 
3.3 The temperature second approximation described in Equation 6 is based on the existence of fonction F that 
translates the detector signal to temperature. 



Equation 7. 



As was explained in paragraph 2.3 the amount of energy exchanged between the detector and the optics is not negligible. 
Assuming that the optics temperature influence can be expressed by a polynomial expansion of order R and assuming 
that the detector signal translation to temperature can be described by a polynomial expansion of order L, the most 
general form of function F is given by: 
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Figure 8. Second order polynomial approximation. 
The normalized error is shown at left, the temperature error at right 
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Figure 9. Third order polynomial approximation. 
The normalized error is shown at left, the temperature error at right 
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Figure 10. Fourth order polynomial approximation. 
The normalized error is shown at left, the temperature error at right 



L R 

F(Sig, To) = (X ^rj • To^)Sig^ [Temperature in degrees Kelvin] 

Equation 8. 

The minimum mean square error algorithm can calculate the constants a^j assuming that the detector's signal has been 

acquired for different view and optics temperatures. The number of linear equations that have to be simultaneously 
solved is (1+ 1/ Although Opgal uses a 32-bit floating point arithmetic unit, the accwacy is not adequate for such 
a task. Therefore we separate the view temperature from the optics influence, histead of having one single F function, 
there are several functions Fro(Sig) that cover the optics ambient temperature range in small steps of 5 degrees. 

FjoiSig) ^j^a^^j ■ Sig' 
Equation 9 

The optics change in temperature is relatively slow. Consequently, the additional work performed by the digital signal 
processor (DSP) in order to calculate, by linear interpolation, the required a^^j constants for the new optics temperature 

is negligible. 

3.5 In order to build a real radiometer, additional variables have to be taken into account, such as the target 
emissivity and ambient temperature influence. However, we would like to mention the influence of the system 
modulation transfer function (MTF) on temperature measurement accuracy. A given object, at a constant temperature, 



will show different temperatures at different distances from the camera. The MTF influence can be compensated, at least 
for some limited spatial frequency range. 




Fourier transform of camera impulse response 

Figure 11. Fourier transform of camera impulse response in cycles per mm, including optics detector and 
electronics. 

Figure 1 1 describes a typical two-dimensional Fourier transform of an uncooled microbolometer camera. In order to 
compensate the MTF degradation versus spatial frequency, an inverse filter to the MTF response was designed and 
software implemented. Figure 12 represents the high pass frequency response of the filter designed in order to 
compensate the MTF degradation. 

3.6 The ideas described in this paper have been implemented on an instrument design to detect human beings suffering 
from fever. The instrument is monitoring the face temperature and the areas in the image above some given temperature 
threshold are flickering. The instrument covers a relatively very small range in temperatures between 30 to 42 degrees 
Celsius within an accuracy of ±0.25 degrees Celsius. 
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Hgh PassFflter Frequency Response. This filter is the MTF inverse 
approxlmatron for a detector that contains 320 by 240 dements. 



Figure 12. MTF inverse filter. 



4 Conclusions 



The basic concept of temperature measurement using a regular uncooled microbolometer detector was described in 
this paper. The additional hardware required for the temperature measurement capability is really limited and 
inexpensive. Opgal expects that the capability to measure the temperature will expand the applications span of the 
regular uncooled microbolometer FLIR cameras. 
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To: 'Geoff Melnick' 
Subject: RE OPG 1 

I really had something like that tn mind, 
claim 3: a flag will be better than the flag 

» 

OK TO FILE IN THE IPO 
in the name of OPGAL Ltd. 

Jonathan Zohn 

I P Manager 



This e-maii message has been sent by EJbit Systems Ltd. or BOp Ltd. and is for the use of the intended 
recipients only. The message may contain infomiation that is client-Forney privileged or confidential. If you 
are not the intended recipient you are hereby notified that any use, distribution or copying of this 
communication is strictfy prohibited, and you are requested to delete the e-mail and any attachments and 
notify the sender immediately. 



— Original i^essage — 

From; Geoff Melnick [mailto:geofF@ipatentco.il] 

Sentt Thursday, August 07, 2003 7:23 PM 

To: Jonathan Zohn [ELOP] 

Subject? RE: OPG 1 

HlJonathan, 

Attached please find my revisions to the claims. 

What I have done is removed the additions from claim 1 and written them out as separate dependent 
claims. That way you get the best of both worlds. It is cleariy written out but doesn't limit. 



Geoffrey. L. Melnick 
Patent Attorney 
EHRLTCH & PARTNERS • 
G. E. Ehrlich (1995) Ltd. 
Ayaion Tower, 15 Ih Floor, 

1 1 Menachem Begin Street, 52 521 Ramat Gan, Israel 

Tel.: 03^6127676 Fax.: 03-61 27575 .Mobile 054962625 Home 08 936 2304 

My address: geoCf@ipacen£.co.il 

Our Website: www, ipatenLco.il 

Alliance: Granot, Strauss, Adar Co. - Website: www.granot-law.co.il 

Tlus e-mail message is for the sole use of the intmKfed recipiem(s) and may contain confldeotial and privileged infonnation. Any 
unauthcnizcd review, use. disdosuce or distribution is prohibited. If you are not the intended recq}ient. please contact the sender by 
reply e-mail and destroy all CjOpies of the original message. 



Thi& e-mail message has been sent by Elop EleotroOptics Industries Ltd. and Is for the use of the 
Intended recipients only. The message may contain mformation that is privileged or confidential.lf 
you are not the intended recipient you are hereby notified that any use, distribution or copying of 
this communication is strictly prohibited, and you are requested to delete the e-mail and any 
attachments and notify the sender Immediately. 
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